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In certain extensions of General Relativity, wormholes generated by spherically symmetric electric 
fields can resolve black hole singularities without necessarily removing curvature divergences. This 
is shown by studying geodesic completeness, the behavior of time-like congruences going through the 
divergent region, and by means of scattering of waves off the wormhole. This provides an example of 
the logical independence between curvature divergences and space-time singularities, concepts very 
often identified with each other in the literature. 
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I. INTRODUCTION 

The blow up of curvature scalars or of components of 
the Riemann tensor is typically used in the literature 
to put forward the existence of space-time singularities. 
Indeed, it is common to use those divergences to argue 
that when curvature scalars get close or surpass Planck- 
ian scales it is time to replace Einstein’s theory by an 
improved description, opening a door to quantum the¬ 
ories of gravity. This intuitive view has even shaped 
numerous approaches in the search of nonsingular space- 
times by t^ing to build theories with bounded curvature 
scalars [lH6|. However, from a formal perspective, the 
most widely accepted criterion to characterize a singular 
space-time is through the existence of incorrmlete paths 
(inextendible time-like or null geodesics) [3-l3 (see also 
[loll and references therein). 

Given that geodesics are classical geometrical enti¬ 
ties associated with idealized (structureless) physical ob¬ 
servers, it is unclear what a quantum theory of gravity 
should say about them. The focus is thus typically turned 
back to curvature divergences as one of the disturbing el¬ 
ements an improved theory of gravity should get rid of. 
A combined approach to better understand the correla¬ 
tion observed between the incompleteness of geodesics 
and the appearance of curvature divergences consists on 
exploring the impact of curvature divergences on physi¬ 
cal observers represented by geodesic congruences. This 
has led to a classification of the strength of curvature 
divergences [IMl according to whether an object go¬ 
ing through them is crushed to zero volume or ripped 
apart by infinite tidal forces (strong case), or if it is not 
severely affected (weak case). In the end, however, one 


must accept that, to the best of our knowledge, matter 
and energy have quantum properties and a fundamental 
wave-like behavior, which demands for a characterization 
of curvature divergences and space-time singularities by 
means of quantum scattering experiments [16 1. 

In this work we provide an example of black hole space- 
time in which a wormhole gets rid of the central singu¬ 
larity without necessarily avoiding curvature divergences. 
This geometry is an exact electrovacuum solution of cer¬ 
tain high-energy extensions of general relativity (GR) for¬ 
mulated in metric-affine spaces (see [T7l - [l^ . where this 
geometry is derived in detail). The requirement of a 
metric-affine geometry is essential to avoid ghosts and 
higher-order derivative equations. The resulting space- 
time roughly consists on two copies of the Reissner- 
Nordstrom or Schwarzschild black hole solution^ con¬ 
nected by a wormhole in a small region near the center. 
This wormhole is supported by an electric field and, as 
such, satisfies all the classical energy conditions. 

Incomplete geodesics that in GR end (or start) at the 
central singularity can now go through the wormhole, 
thus avoiding incompleteness (20|. The wormhole throat, 
however, generically exhibits divergent curvature scalars. 
Nonetheless, we find that physical observers (described as 
time-like congruences) do not experience a ny pathologi¬ 
cal behavior upon crossing the wormhole [^. We also 
find that a scalar held propagating in this background is 
everywhere well-behaved, and compute the transmission 
coefficients and cross section of scalar waves as they in¬ 
teract with the wormhole [l^. The resolution of black 
hole singularities in this model provides new insights on 
how this disturbing aspect of classical gravitation may 
be cured at high energies. 
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^ Depending on particular parameters, two copies of Minkowski 
are also possible. 
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II. BACKGROUND GEOMETRY AND 
GEODESICS 


The geometry we are going to study arises when a 
static, spherically symmetric electric field is coupled to a 
certain extension of GR (see [Hill for details). In fact, 
this solution appears in two different gravity theories, 
namely, the quadratic gravity model with Lagrangian 
i? — I + aR^) and the Born-Infeld gravity model 
“ '/-^)- The line element, whose 
derivation was worked out in mill , takes the form 

ds^ = —A{r)dt^ H- \-^dx^ + r^{x)dQ^ , (1) 

B{r) 


with dVi^ = dO"^ + sin0^d(/3^, B{r) = A{r)( 7 ^, and 


’■{x) = 


+ yjx'^ + 4r4 


A{x) = — 
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1 - 


rs (1 + <5iG(r-)) 
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( 2 ) 

( 3 ) 

( 4 ) 


where a; s] — oo,+oo[, but r{x) > Tc > 0. The function 
G(z), with 2 ; = r/rc, can be written as 


G{.z) = 1 [fs/iiz] + fr/iiz)] , (5) 


where fx{z) = 2 T’i[^,A, |,1 — z'^] is a hypergeometric 
function, and 6 c ~ 0.572069 is a constant. By careful in¬ 
spection, one can verify that this line element describes a 
wormhole [ 2 ^. To see this, Eq.Q is particularly useful, 
as it shows that the area function A = Airr'^ {x) has a min¬ 
imum Amin = 47rr^ at a; = 0, where the wormhole throat 
is located. Here Xc = 4 is a length scale character¬ 

izing the departures from GR HP (in fact e = —21^), 
= 2Gq^ is the charge radius, rg the Schwarzschild 

radius, and di = { 2 rs)~^ is the charge-to-mass 

ratio, which plays an important role in the classification 
of the solutions. 

As already noticed, the geometry described above is 
an exact solution of certain extensions of GR such as the 
Born-Infeld model and quadratic gravity [13 coupled 
to an electric field. The reason why such solutions were 
not known previously is because we considered those the¬ 
ories from a metric-affine perspective, i.e., assuming that 
the metric and the affine connection are independent geo¬ 
metric entities, which yields field equations different from 
those found in the typical metric approach This 

formulation of modified gravity yields second-order field 
equations that recover Minkowski (or de Sitter) space- 
time in vacuum, yields no extra propagating degrees of 
freedom, and is ghost free. The matter source threading 
the corresponding space-time is an electrostatic spher¬ 
ically symmetric Maxwell field, which means that the 
classical energy conditions are all satisfied. We also note 


that this static geometry can be obtained by considering 
the dynamical process of collapse of charged radiation 
fluids in a Vaidya-type space-time [^ . 

From the above definitions, one can verify that the 
GR geometry is quickly recovered as soon as |a;| is 
slightly greater than Xc (recovering B{r) Ri A{r) Ri 
1 — rs/r + rg/{2r^), and ~ x^) and also in the 
limit Ic —>■ 0. Thus, for scales larger than the worm- 
hole size Tc, this space-time looks much the same as the 
standard Reissner-Nordstrom solution of the Einstein- 
Maxwell field equations. As x —>■ 0, however, the area 
function r^(x) ~ -|-x^/2 reaches a minimum, and A(r) 

dramatically departs from GR, behaving as (from now 
on we set le = Ipianck for notational convenience) 


A{r) 


1 _ IS I _:_2_ 
r ' 2r2 

ANc SiSc 



Na-Nq 

2Nc 


if |x| » Tc 

if X —0 


( 6 ) 


where Nq = \q/e\ is the number of charges and Nc ~ 
16.55. This behaviour turns the Kretschmann scalar from 
R^^is^xRa^^^ ~ 14r^/r® in the Reissner-Nordstrom solu¬ 
tion into R°‘pjxRa^'^^ ~ (<^1 - 6 c)^K 2 l{r - r^)® -b ((5i - 
Sc)Ki/{r — rc)^G^-XQ in our case, with the Ki constants. 

If the charge-to-mass ratio i5i takes the particular value 
i5i = Sc, the metric and all curvature scalars are regular 
everywhere and, naturally, all time-like and null geodesics 
smoothly extend across the wormhole. In this case, if 
Nq > Nc, which implies lim 2 ;^o Air) < 0, the worm- 
hole is hidden in between two event horizons symmet¬ 
rically located around x = 0 (r = Tc), connecting in 
this way two Schwarzschild-like geometries. If Nq < Nc, 
there are no event horizons, and the wormhole connects 
two Minkowski-like spaces. When (5i y/ Sc, however, cur¬ 
vature divergences arise at the throat and the geodesic 
equation must be studied in detail to determine if those 
configurations are geodesically complete or not. 

In terms of an affine parameter A, the geodesics of 0 
satisfy [ 2 ^ (k = 0,1 for null and time-like, respectively) 


1 





A{r) 



( 7 ) 


where E and L are constants which, in the time-like case, 
represent the total energy and angular momentum per 
unit mass. For null geodesics with L = 0, we get 



whose solution is plotted and compared with the GR pre¬ 
diction in FiglTJ Given that ([8]) is independent of Si, all 
configurations (with or without curvature divergences) 
have an identical behavior, which confirms the extendibil- 
ity of null geodesics across x = 0 in these space-times. 



















3 


A 



Figure 1. Affine parameter A(a:) as a function of the radial 
coordinate x for radial null geodesics (outgoing in a; > 0). In 
the GR case (green dashed curve in the upper right quadrant), 
A = X is only dehned for x > 0. In this plot E = 1 and the 
horizontal axis is measured in units of r^- 


For null and time-like geodesics with L ^ 0, near the 
wormhole the effective potential 


V{x) = A 



(9) 


can be approximated as V(x) —a/\x\, with a = 

readily sees that if 5i > 6c 
(Reissner-Nordstrom-like case), an infinite potential bar¬ 
rier arises which prevents any geodesic from reaching 
X = 0. This behavior is also observed in GR, where such 
geodesics can never reach the singularity (26|. In the 
Schwarzschild-like case, (5i < 5c, the potential becomes 
infinitely attractive as x —>■ 0 and turns © into 


dx 


±1 

2 


A(x) 




( 10 ) 


where ± is for outgoing/ingoing geodesics. In the case 
of GR, time-like radial geodesics near the singularity be¬ 
have as A(r) Ri Given that r > 0, ingo¬ 

ing/outgoing geodesics in GR end/start at A(0) = 0 and, 
therefore, are incomplete. In our case, on the contrary, 
the fact that x g] — oo, -|-oo[ allows to smoothly extend 
the affine parameter A across x = 0 to the whole real axis. 
It is also worth noting that the conclusion of the singular¬ 
ity theorems [l^, that in GR prevent the Schwarzschild 
geometry to be extended, does not hold for these worm- 
hole space-times, because although the electromagnetic 
stress-energy tensor fulfills the weak energy condition, 
> 0 for k^^ non space-like, it does no longer 
imply R^,j{g)k^^k'^ > 0 in our gravity theories. 


III. CONGRUENCES 

We have just seen that there are no restrictions for null, 
time-like or space-like geodesics to be extended across the 


wormhole (when they can reach it). Nonetheless, given 
that physical observers can be more accurately described 
as geodesic congruences, one should consider the im¬ 
pact of curvature divergences in the evolution of nearby 
geodesics. For this purpose, we consider the line element 
© written in coordinates adapted to an observer in free 
fall and study if causal contact between nearby time-like 
observers is lost at any time of the transit across the hole. 
In freely falling coordinates, o turns into [2l| 

ds^ = -dX^ + -f r^{\, C}dXf , (II) 

where A is the proper time (affine parameter), ^ mea¬ 
sures the radial separation between geodesics, and = 
dy/dX, with dy = dx/{l + r‘^/r^{x)). Focusing on the 
Schwarzschild-like case, < Sc, which is the only config¬ 
uration in which time-like geodesics can go through the 
wormhole, (m^)^ can be approximated as (m^)^ ~ — 

(||A — E^\)~i, which turns (fTTl) into 




/3 


( 12 ) 


One should note that as the wormhole throat is ap¬ 
proached, X — ^ 0, the physical spatial distance be¬ 

tween any two infinitesimally nearby geodesics diverges: 
dlphys = (f|A —(see [2lj for a more elab¬ 
orate discussion of this point using Jacobi fields). This 
could be seen as an indication that strong tidal forces 
could rip apart any infalling body, thus signaling the 
presence of a strong singularity. However, for any finite 
comoving separation the physical spatial 

distance Iphys = f is always finite. 
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Phys 


/a> 

1/3 1 
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|A-i?^o|2/3-|A-£;a|2/3 


(13) 


which casts doubts on the fate of such bodies and the 
actual strength of those curvature divergences. 

It is thus necessary to clarify if the constituents mak¬ 
ing up an object that reaches the wormhole lose causal 
contact because of the divergent spatial stretching that 
affects its infinitesimal elements (infinite tidal forces). In 
such a case, the interactions that keep the object cohe- 
sioned would no longer be effective, resulting in disrup¬ 
tion or disintegration of the body. However, as shown in 
Fig 121 a fiducial observer at ^ = 0 never loses causal con¬ 
tact with its neighbors. Indeed, in [2l| we show that the 
proper time a light ray takes in a round trip from ^ = 0 
to any nearby geodesic is always finite as the wormhole is 
crossed. The divergent stretching of (infinitesimal) spa¬ 
tial distances is, therefore, not observable. We can thus 
conclude that extended objects going through the worm- 
hole, where curvature scalars generically diverge, do not 
undergo destructive deformations and, therefore, can ef¬ 
fectively cross this apparently troublesome region. 
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Figure 2. Trajectories of light rays emitted by a freely falling 
observer from ^ = 0 at different times shortly before reaching 
the wormhole throat in the Schwarzschild-like conhguration. 
The rays going to the left/right represent ingoing/outgoing 
null geodesics. Given that the observer is inside an event 
horizon, both ingoing and outgoing light rays end up hitting 
the wormhole, which is located at the oblique (solid black) 
line A — = 0 (in the plot E = 1, a = 3). 


IV. SCATTERING EXPERIMENTS 


As a third test to verify that curvature divergences 
do not affect the well-posedness of physical laws in our 
space-time o, we consider the propagation of scalar 
waves in this background. According to the cosmic cen¬ 
sorship conjectnre, naked singularities should not exist 
in Nature. This idea is traditionally invoked in the lit¬ 
erature to interpret naked singularities as unphysical ar¬ 
tifacts of Einstein’s theory. From this physical perspec¬ 
tive, therefore, one can say that naked divergences are 
the worst-case scenario. From a technical viewpoint, this 
turns out to be the simplest situation because there ex¬ 
ists a time-like Killing vector over the whole space, which 
facilitates the choice of coordinates and the separation of 
variables. Our approach is also valid for the region inside 
the inner horizon, where the wormhole is a time-like hy¬ 
persurface. We will thus focus on configurations with 
i5i > Sc (Reissner-Nordstrom-like). Taking the scalar 
field equation (□ — m?)(j) = 0, we decompose it in modes 
of the form (j)uj,ira = e~'^‘^*Yim{d,<p)fu:,i{x)/r{x), where 
Yim{d,^) represents spherical harmonics. Using the ra¬ 
dial coordinate y = f dx/A{l 4- the fui,i{x) are 

governed by a Schrodinger-like equation of the form 


where 


fvv + ^efff — w /, 

+ A{r) 


( 14 ) 


(15) 


This potential quickly converges to the GR prediction for 
r ^ Tc but near r = Vc diverges as 14// ~ k/\y\^^^, where 


k = 


' (5i - 6c)N, {Nc[m\l + 1 + 1(1 + 1)] - N,) 


SiScNc 


Ac(8r3)i/2 


(16) 


Low frequency modes are almost entirely reflected by 
the centrifugal barrier in much the same way as it hap¬ 
pens in GR (see [l^ for details). We thus focus on high- 
energy modes with sufficient energy to overcome this bar¬ 
rier and interact with the wormhole. The effective po¬ 
tential is divergent at the wormhole throat, but as is well 
known from the example of the Coulomb potential, this 
does not mean that the wave function is singular. Actu¬ 
ally, the leading behaviour of the wave function ^ y^ 
is given by the characteristic exponents p = 0,1, which 
are perfectly regular. Considering the m? = 0 case for 
simplicity, it is easy to see that 14// may transit from 
an infinite well to an infinite barrier as I increases if 
Nq > Nc- The divergent barrier is reminiscent of the 
repulsive potential seen by geodesics in this same config¬ 
urations. Using the mode decomposition, we thus define 
an incoming wave packet from past null infinity (in the 
naked case, or from the inner horizon in the black hole 
case) and study its interaction with the wormhole. As 
noticed above, the behavior will depend dramatically on 
the angular momentum of the incident mode. 

With a notational change, the information contained 
in the parameters l,Nq,di and the frequency lu can be 
condensed in a (dimensionless) parameter a dehned as 

y' = \k\3y a = \k\~iuj . (17) 

The wave equation in the relevant region then becomes 

fy,y, + ± ^)/ = 0 , (18) 

with the ± sign set by the sign of k (infinite well or 
barrier). Now we can compute numerically the trans¬ 
mission coefficient for a given value of a. The full po¬ 
tential decays appropiately as the radius increases and 
the wave function can be approximated asymptotically 
as Ain exp(ia;(t — r)) + Aout exp(iuj(t + r)) on one side of 
the wormhole, and Bin exp(zaj(i — 'r))-|-i?out exp{iu!{t+r)) 
on the other. To obtain the transmission coefficient, we 
give initial conditions /, /' at a point Xout across the 
wormhole with the condition Bin = 0, and evolve the 
wave function up to a point Xin on the side of the worm- 
hole where the wave was sent. From the values of f(xin), 
f'(xin), f(xout), f'(xout) we can extract the amplitudes 
Ain, Aout, Bout- The transmission coefficient is simply 
T = I Bout/Ain p. There is a a caveat if we use the ap¬ 
proximated effective potential, because it is long-range 
and the wave function cannot be approximated with ex¬ 
ponentials as before. We can use the WKB approxima¬ 
tion far from the wormhole to obtain the amplitudes Ain, 
Aout, Bout, but this would introduce a phase respect to 
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Figure 4. Transmission cross section versus frequency for a 
naked divergence with = 1, Nq = 10, calculated numeri¬ 
cally (dots) for the exact 14//. The continuous line shows the 
approximation a cx uj 2 . 


the amplitudes calculated using the full potential. How¬ 
ever, the transmission coefficient would not be affected, 
and this is the method we have used. 

In Fig. [3] the transmission coefficient as a function of 
a is plotted (for the approximated 14//). For k < 0 
(infinite well), most of the wave is transmitted, with the 
transmission coefficient tending to 1 as a grows. For 
k > 0 (infinite barrier), a sigmoid profile arises like in 
typical barrier experiments. For a above ^ 1.5, the wave 
is transmitted almost entirely, and below that threshold 
it is almost completely reflected. 

Consider now a case with constant w. Depending on 
the number of charges, k can be either positive or neg¬ 
ative for I = 0. With growing I, k grows as ~ P. The 
transmission for k negative is very high, and approaches 
1 as fc —>■ 0. When k changes sign as I grows, a —>■ 00 and 
the transmission is near 1. For bigger Z’s, a begins to de¬ 
crease, and around the threshold a ~ 1.5, characterized 
by I = 4iaxj the partial waves change from being almost 
totally transmitted to being almost totally reflected. A 
rough estimate of the cross section can thus be obtained 
by considering that the transmission is 1 for partial waves 
with I < ^inax, and 0 for I > /max- The transmission cross 
section would thus be 


CT — —^ {21 -|- 1)1 — ^(1 + ^max) (19) 

1=0 

The variation of a in terms of w is illustrated in Fig. |4l 
As w —>■ 00 , we have Zmax oc which confirms that a 
decreases as ~ a;“ 2 . 


V. SUMMARY AND CONCLUSIONS 

Two fundamental properties typically demanded from 
quantum theories of gravity (often seen as one and the 
same) are the removal of space-time singularities and 
the regularization of classical curvature divergences. We 
have shown here that only the former is physically rele¬ 
vant and can be achieved in classical geometric scenar¬ 
ios with nontrivial topology and satisfying all the energy 
conditions. Though our analysis has focused on a four¬ 
dimensional setting, our results are still valid in arbitrary 
dimensions, where exact analytical solutions with worm- 
hole structure, and which are geodesically complete, have 
been found (we note that geodesic completeness via 
wormholes can also be achieved in certain metric-affine 
/(i?) theories [1^). Higher dimensional models confirm 
that the resolution of black hole singularities is a highly 
non-perturbative phenomenon, being the GR solution an 
excellent approximation down to the very throat of the 
wormhole, where the geometry quickly changes to ac¬ 
count for this topological structure. We also note that 
the theories that generate the solution o also replace 
the big bang singularity by a cosmic bounce in a very 
robust manner |3fll432j | . 

The generic absence of space-time singularities in these 
theories without requiring exotic matter-energy sources 
puts forward the existence of an important gap in our 
understanding of gravitation in metric-affine geometries. 
Interestingly, this type of geometries seem to play an 
important role in the description of continuum systems 
with a microstructure [ 3 ^ . such as Bravais crystals or 
graphene. This provides further motivation for their 
study in gravitational scenarios. 
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